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NOTATION 


Symbol  Definition  Dimensions 

2 

A  Submerged  cross-sectional  area  ft 

Ag,  A^  . . .  Coefficients  in  ship  line  function 

a^,  a^,  a,.,  m,  Parameters  used  in  determining  added 

_  mass;  see  Appendix  B2.3 

n,  p,  ft ,  r  ,  r2, 

sr  S2'  f 


dx  _  d£  Longitudinal  velocity  ft /sec 

dt  dt 

E  Modulus  of  elasticity  in  tension  and 

compression  ton-ft^ 

El  Flexural  rigidity  of  hull  ton-ft^ 

f  (y)  Half -breadth  corresponding  to 

level  y  ft 

2 

G  Modulus  of  elasticity  in  shear  ton/ft 

g  Vertical  acceleration  of  hull  due  to  2 

gravity  ft /sec 

h  Wave  height  ft 


v 


Symbol 


Definition 


Dimensions 


Moment  of  inertia  of  submerged  cross 
section  about  transverse  axis  in  free  ^ 

surface  ft 

Function  of  I  defined  in  Appendix  B2.3  dimensionless 


Mass  polar  moment  of  inertia  per  unit 
length  of  section  of  hull  A  x  long  taken 
about  a  horizontal  axis  through  its  c.g. 
(including  allowance  for  virtual  mass 
of  surrounding  water) 

Defined  in  Section  B  2.  6 
Shear  flexibility  factor; 


r(,-f) 


Cos  ft 


ton-sec 


ICAG 

Shear  rigidity  of  hull 

ton 

M 

Bending  moment  on  hull 

ft  -ton 

m 

Hydrodynamic  added  mass  per  unit 
length 

slugs /ft  = 

iu  2/ -*2 

lb-sec  /it  ; 
also 

.  2 
ton-sec 

ft2 

m 

s 

Ship  mass  per  unit  length 

slugs  / ft;  also 
ton-sec  /ft 

n,  n  +  1/2 

Station  and  midstation  numbers, 
respectively,  for  ship  where 
stations  run  0  (stern),  1,  ...  N-l, 

N  (bow)  and  midstations  run  1/2, 

1  1/2,  ...  N  -  1/2;  N  is  also  distance 
along  normal  to  surface 

ft 

P(x,  t) 

Total  force  per  unit  length  acting  upon 
the  ship  hull 

ton/ft 

P  ,  P 
e  e 

Defined  on  Pages  70-71  of  Reference  1 

VI 
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Symbol 

Definition 

Dimensions 

t 

Time 

sec 

U 

Uniform  forward  velocity  of  ship 
along  its  heading 

ft/sec 

u 

In  general,  the  horizontal  component 
of  particle  velocity  at  any  depth  in  the 
undisturbed  fluid 

ft/sec 

V 

Vertical  shearing  force  acting  on  hull 

ton 

V 

w 

Velocity  of  wave  propagation 

ft/sec 

V,  V 

In  general,  the  vertical  component  of 
particle  velocity  and  acceleration, 
respectively,  at  any  depth 

ft /sec  and 

ft/sec^, 

respectively 

W 

Weight  of  hull  per  unit  length 

ton /ft 

X 

Abscissa  or  axis  of  abscissas  in  a 
rectangular  coordinate  system;  a 
coordinate  indicating  distance  from 
wave  crest  of  point  on  ship;  horizontal 
displacement  of  particle  at  time  t; 
particle  coordinate  of  particle  at  rest; 
longitudinal  distance  measured  from 
stern  of  ship 

ft 

X 

C.  g. 

Distance  from  stern  to  center  of 
gravity 

ft 

(y,  z) 

Coordinate  of  ship  profile  for  each 
section 

ft 

y 

Ordinate  or  axis  of  ordinate  in  a 
rectangular  coordinate  system; 
immersion;  particle  coordinate  of  a 
particle  at  rest 

ft 

y 

Lateral  deflection  of  hull 

ft 

y 

Immersion  velocity 

ft/sec 

Vll 


Symbol 

Definitions 

Dimensions 

y*.  y",  y'" 

Draft;  (See  Appendix  B  of  Reference  1) 

ft 

Elastic  (flexural)  component  of  y 

ft 

V  K  - 

’Displacement  and  vertical  velocity  of 
keel,  respectively,  due  to  heave 

ft  and  ft /sec, 
respectively 

yi 

i  +  1  water  level,  i  =  0, 

1  ...  N  -  1 

ft 

V 

Displacement  and  vertical  velocity  of 
keel,  respectively,  due  to  pitch 

ft  and  ft /sec, 
respectively 

yr 

Downward  relative  velocity  of  hull  with 
respect  to  fluid  undisturbed  by  ship 
(not  wave) 

ft /sec 

yw 

Elevation  of  wave  surface  above  mean 
water  level 

ft 

Z,  z 

Complex  variables  for  the  transformed 
and  circle  planes,  respectively 

0 

Deadrise  angle 

deg 

y>  y 

Rotation  and  angular  velocity,  respectively, 
of  transverse  sections  with  respect  to  a  radians  and 
horizontal  axis  radians /sec, 

respectively 

At 

Time  interval  or  increment 

sec 

A*  x 

Length  of  element,  Jtl 20;  increment 
of  length 

ft 

c 

Distance  of  station  n  forward  of  position 
of  heave  meter 

ft 

• 

V>  7 

Vertical  displacement  from  rest  condition 
and  associated  vertical  velocity  of 
particle  at  arbitrary  depth  y;  section  area 
coefficient  (dimensionless) 

ft  and 

ft/sec, 

respectively 

vau 


Symbol 


Definition 


Dimensions 


9  2  n  x 

A 

Q  Value  of  0  at  any  station  n  (equal  to 

0.1366  ■,  ■-1~AA  g  +  0.  2732  n-  2.732) 
0.1096 

q  Ship's  heading 

s 

*  Wave  length 

H  Mass  per  unit  length  of  hull  (including 

allowance  for  virtual  mass) 

• 

£  is  distance  of  any  station  n  along 
keel  forward  of  midlength;  horizontal 
velocity  of  a  particle  at  arbitrary  depth 

y 

>  Value  of  £  at  station  n 

€  n 

p  Density  of  sea  water 

^ » 'P  Pitch  angle  and  angular  velocity, 

respectively 

w  Angular  velocity 


radians 

radians 

deg 

ft 

ton-sec^/ft2 

ft  and  ft  /sec, 
respectively 

ft 

2  4 

lb -sec  /ft 

radians  or  deg 
and  deg /sec, 
respectively 

radians /sec 
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ABSTRACT 


This  report  describes  a  method  for  obtaining  digital 
computer  solutions  for  the  excitation  forces  on  and  transient 
response  of  a  ship  subject  to  (hydrodynamic)  slam  when 
certain  basic  data  are  obtained  by  computation  rather  than 
by  measurement. 

The  method  is  based  on  a  theory  shown  in  TMB 
Report  1511  to  be  in  good  agreement  with  experiment  for  a 
particular  application. 

INTRODUCTION 

A  theoretical  analysis  and  manual  computation  of  the  "slamming" 
(hydrodynamic)  forces  acting  on  a  ship,  based  upon  an  experimental 
knowledge  of  the  rigid -body  motions  of  the  ship  relative  to  the  wave,  was 
presented  in  Reference  1.  These  forces  are  considered  to  be  due  to 
the  time  rate  of  change  of  fluid  momentum  and  to  buoyancy  forces  incident 
to  immersion  of  the  hull.  In  addition,  by  use  of  a  digital  computer,  a 
computation  was  made  of  the  transient  elastic  response  and  associated 
hull  girder  stresses  of  the  ship  due  to  the  total  force  exerted  by  the  fluid 

In  this  report  we  are  considering  the  slamming  forces  associated  with 
immersion  of  a  flared  bow  rather  than  bottom  impact  or  side  pounding 
forces.  The  forces  referred  to  here  are  the  time-varying  equivalent  or 
integrated  load  force  on  each  cross  section  rather  than  the  detailed  load 
distribution  around  the  cross  section. 

References  are  listed  on  page  37. 


on  the  ship.  A  comparison  between  the  theoretical  and  measured  stresses 
for  the  Dutch  destroyer  showed  good  agreement.  This  method,  therefore, 
offers  promise  for  evaluating  various  hull  shapes  in  the  early  design 
stage  if  a  practical  computational  procedure,  requiring  a  minimum  of  in¬ 
put  data,  can  be  established.  Specifically,  it  is  desirable  to  reduce  the 
complexity,  expense,  time,  labor,  and  errors  often  involved  in  making 
numerical  computations  for  predicting  the  slamming  forces  on  and  the 
dynamic  (elastic  and  rigid -body)  response  of  a  ship  girder  without  a 
prior  experimental  knowledge  of  the  rigid -body  and  flexural  motions  of  a 
ship  relative  to  the  wave. 

Accordingly,  the  objective  of  this  report  is  to  describe  a  method 
for  obtaining  digital  computer  solutions  for  the  excitation  forces  on  and 
transient  response  of  a  slammed  ship  when  certain  basic  data  are  obtained 
by  computation  rather  than  by  measurement.  Moreover,  as  shown  in 
Reference  2,  those  parts  of  the  basic  data  which  involve  ship  properties 
may  also  be  computed  digitally  rather  than  manually  using  tabulated  quan¬ 
tities  systematically  obtained  in  a.  prescribed  fashion  directly  from  ship 
plans. 

A  digital  computer  (IBM  7090)  is  presently  being  coded  at  the 
Applied  Mathematics  Laboratory  of  the  David  Taylor  Model  Basin  to  ob¬ 
tain  solutions  in  accordance  with  this  method;  this  code  includes  a  digital 
computer  routine  for  computing  the  hull  parameters. 
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METHOD  OF  ATTACK 


The  following  procedure,  based  upon  the  theory  of  Reference  1,  is 
used  to  solve  for  the  excitation  forces  on  and  the  response  of  the  ship  for 
an  arbitrary  transient  wave  of  encounter: 

a.  The  basic  finite  difference  equations  of  motion  used  to  obtain 
the  damped  transient  vertical  response  of  a  ship  to  slamming  forces, 
derived  in  Appendix  F  of  Reference  1,  are  coded  for  solution.  For  con¬ 
venience  of  reference,  the  corresponding  differential  equations  are  re¬ 
written  in  Appendix  A. 

b.  Time -independent  ship  data  required  for  solution  of  these 

equations  are  furnished  to  the  digital  computer;  see  Appendix  B. 

c|x  *  __  ^ 

c.  Time -dependent  data  (e.g.,  m,  ,  yr  ,Pg,  c,  etc.)  required 

for  solution  of  these  equations,  are  computed  by  the  digital  computer  from 

other  time -independent  data,  some  additional  constants,  and  several 

3 

points  that  describe  each  ship  line  (i.e.,  cross-sectional  profile);  see 
Appendix  B. 

d.  The  dynamic  elastic  and  rigid -body  response  of  the  ship  due 
to  slam  is  represented  by  the  digital  computer  solution  of  the  finite  dif¬ 
ference  equations  corresponding  to  the  differential  equations  of  Appendix 
A,  using  data  obtained  in  accordance  with  Items  b  and  c. 

An  alternative  method  of  attack  which  was  presented  in  Reference 
1,  is  compare/!  with  this  new  method  in  Appendix  C. 
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The  application  of  this  method  to  an  arbitrary  transient  wave  of 
encounter  is  discussed  in  Appendix  D. 

DISCUSSION 

In  the  process  of  digitally  computing  the  forces  on  and  transient 

response  of  a  slammed  ship,  intermediate  calculations  of  added  (i.e., 

virtual)  mass  and  ship  properties  are  performed  by  means  of  a  digital 

computer.  The  digital  computation  of  the  added  mass,  using  functions 

3 

approximating  a  ship  line  (i.e.,  here  a  cross-sectional  profile)  and  of 
4  5 

Landweber's  J  equations  is  given  in  Appendix  B.  Landweber's  equa- 

6  7 

tions  include  those  of  Lewis  and  Prohaska  as  a  special  case.  Refer¬ 
ence  2  describes  the  digital  computation  of  the  ship  properties  (i.  e.,  hull 
parameters)  using  data  tabulations  obtained  by  a  pre-established  orderly 
procedure  from  ship  plans.  These  apparently  novel  methods  for  obtaining 
virtual  mass  and  ship  properties  should  tend  to  reduce  the  cost,  time, 
complexity,  labor,  and  errors  usually  involved  in  making  such  computa¬ 
tions. 

The  digital  calculations  are  of  special  interest  because  their  ap¬ 
plicability  extends  beyond  the  present  problem  to  other  vibration  problem 
areas  in  which  the  hull  is  also  treated  as  a  beam.  Thus  their  application 

g 

pertains  to  the  determination  of  the  normal  mode  response  of  a  hull, 

9 

steady -state  response  of  a  hull  to  a  sinusoidal  driving  force,  flutter 


response  of  a  hull -control  surface  system  subject  to  hydrodynamic 
forces  on  the  rudder,  ^  etc.  Such  mechanization  fits  the  trend  to¬ 
ward  routinizing  complex  calculations  which  lead  to  eventual  design 
utility. 

CONCLUSIONS 

A  method  has  been  devised  for  obtaining  digital  computer  solutions 
for  excitation  forces  on  and  transient  response  of  a  ship  subject  to  hy¬ 
drodynamic  "slamming"  forces  when  certain  basic  data  are  obtained  by 
computation  rather  than  by  measurement  .  This  method  reduces  the 
complexity,  expense,  time,  labor,  and  errors  heretofore  connected  with 
obtaining  such  solutions.  Consequently,  while  a  ship  is  still  on  the  draw¬ 
ing  board,  it  now  appears  possible  and  practical  to  use  a  high-speed  and 
flexible  computer  to  determine  the  response  to  slam  for  several  different 
combinations  of  parameters  that  represent  different  ship -sea  character¬ 
istics  (see  Appendix  D)  and  operating  conditions. 


See  first  footnote  on' page  1. 
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RECOMMENDATIONS 


It  is  recommended  that  this  method  be  used  in  the  early  design 
stages  to  evaluate  the  response  of  various  hull  shapes  to  slamming 
forces  of  the  type  described  in  this  report.  Computer  results  obtained  by 
systematically  varying  sea  states,  ship  parameters,  ship  geometry  (hull 
shape),  speed,  heading,  etc.,  presented  in  the  form  of  design  curves, 
should  prove  helpful  in  furnishing  useful  hull  design  information  for  direct 
use  by  the  Preliminary  Design  Branch  of  the  Bureau  of  Ships. 
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APPENDIX  A 


DIFFERENTIAL  EQUATIONS  FOR  OBTAINING  RESPONSE 
OF  A  SHIP  TO  SLAMMING  FORCE 


The  basic  differential  equations  used  to  describe  the  motions  of  a 
ship  subject  to  slam,  derived  in  Appendix  F  of  Reference  1,  are  (see 
notation):  (The  corresponding  finite  difference  equations'*'  have  been  coded 
for  the  IBM  7090). 


V-  (x,  l) 


dy(  (x,t) 

Tt 


+  /c(x,t)y(  (x,  t)  +m(x,t) 


dx(x,t) 

dt 


o  - 


1 


KAG  (x) 


dV(x, 
dt  -» 


dV(x,t)  ,  dx(x,  t)  , 

+  — r — —  =  P(x,  t)  +  m(x,t)  - - -  i f)(l) 


dx 


dt 


1  dM(x,t)  d-y(x,t) 


El(x)  dt 


dx 


=  0 


(x) 


dy  (x,t) 
dt 


dM  (x,t) 
dx 


+  V(x,t)  =  0 


l  d  V (x,t)  d%(x>0 

KAG  dt  +  dx 


-  Y(x,t)  =  -ip  (t) 


dyf(x,t) 

dt 


y(  (x,t) 
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where 


p(x,t)  =  m  (x,t)  +  ms  (x) 


c(x,t)  =  B  (x,t)  + 


dm(x,t) 


dm  (x,t) 
+  ~Tx 


dx  (x,t) 

Tt 


P  (x,t)  =  P  (x,t)  -  m 


<c- 


yP  (x,t) 


at 


Pe  (x,t)  =  Pe  (x,t)  -  gms  (x) 
d 

p  (x,t)=——  (my  )  +  (g  +  v)  pA 
e  at 


and  where  for  present  and  future  use  we  define: 
t  Time;  sec 

x  Distance  coordinate  along  longitudinal  centerline  axis  of  ship 

hull;  the  independent  variable  x  as  used  in  the  equations  lies 

9 

along  the  same  axis  as  f  ,  hence  dx  =  d£  and  _  —  ,  ,  . 

dx  d  ^ 

Note,  however,  that  this  is  not  the  same  x  used  in  Figure  1  of 
Reference  1  and  in  other  parts  of  that  reference  where  x  is 
normal  to  the  wavefront;  ft 

The  damping  terms  c  and  c  here  replace  c  and  c  in  the  corresponding 
equations  in  Appendix  FI  of  Reference  1.  In  contrast  to  Reference  1,  this 
change  in  notation  now  permits  a  unique  definition  for  the  symbol  c  ;  i.e., 
half -breadth  of  the  ship  at  a  particular  immersion. 

For  more  complete  definitions,  see  Reference  8. 
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Total  vertical  displacement  of  point  of  beam  (hull)  initially  on 
x-axis,  i.e.,  displacement  of  point  from  equilibrium  position  in 
still  water;  ft 

Elastic  (flexural)  component  of  y;  ft 

Displacement  and  vertical  velocity  of  keel,  respectively,  due  to 
heave;  ft  and  ft/sec,  respectively 

Displacement  and  vertical  velocity  of  keel,  respectively,  due 
to  pitch;  ft  and  ft /sec,  respectively 

Downward  relative  velocity  of  hull  with  respect  to  fluid  undis¬ 
turbed  by  ship  (not  wave);  ft/ sec 

Rotation  of  transverse  section  about  an  axis  through  the  section 
normal  to  the  x-y  plane;  radians 
Bending  moment;  ft -ton 

2  2 

Hydrodynamic  added  mass  per  unit  length;  ton-sec  /ft 

2  2 

Ship  mass  per  unit  length;  ton-sec  /ft 
Net  shear  force  in  y-direction;  ton 

Apparent  mass  per  unit  length  including  allowance  for  virtual 

2  2 

mass  of  surrounding  water;  ton -sec  /ft 

2 

Damping  force  per  unit  velocity  per  unit  length;'  ton-sec /ft 

Total  force  per  unit  length  acting  upon  ship  hull;  ton /ft 

Mass  moment  of  inertia  of  hull  per  unit  length  with  respect  to 

2 

an  axis  normal  to  the  x-y  plane;  ton-sec 

2 

Bending  rigidity;  ton-ft 


KAG  Shear  rigidity;  ton 

2  4 

p  Density  of  sea  water;  lb -sec  /ft 

2 

A  Submerged  cross  sectional  area;  ft 

W  =  m  g  Gravity  load  per  unit  length  where  m  is  the  ship  mass  per 
s  s 

unit  length;  ton /ft 

tit  Pitch  angle  and  angular  velocity,  respectively;  radians  and 

radians /sec 

v,  v  In  general,  the  vertical  component  of  particle  velocity  and 

2 

acceleration,  respectively,  ai  any  depth;  ft/sec  and  ft/sec  , 
respectively 

In  Equations  [Al]  through  [A 5]  the  following  restrictions  arc 

imposed: 

■ 

1  >o 

KAG  (x) 

H  (x,t)  >0 
1 

- >0 

El(x) 

The  ship  is  assumed  to  have  free  ends.  Therefore,  the  end 
(boundary)  conditions  imposed  for  all  times  are: 

V(0,t)mV(l  (0,t)  =  '/  =  0 
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This  is  equivalent  to  V  =■—  =0  at*  =0  ,  and*  =x..for  all  times. 

dx 

The  initial  conditions  required  are: 

y(  (x,o);  y(  (x,o);  y  (x,o);  .1/  (x,o);  V  (x,o) 

The  initial  conditions  actually  supplied  are: 
y(  (x,o);  y(  (x,o);  y  (x,o)  y  (x,o) 

Then  the  required  1/  (x,o)  and  V  (x,o)  are  obtained  by  using: 

By  (x,o) 

M  (x,o)  =  E1  (x)  - 

dx 

d.\l  (x,o) 

V  (x,o)  =  - 

B  x 

In  addition  to  these  end  conditions  and  initial  conditions,  certain 
time -independent  and  time -dependent  data  are  to  be  furnished  as  input  to 
the  computer  (see  Appendix  B). 

The  method  of  computation  follows  the  general  procedure  given  on 
pages  70-71  of  Appendix  FI  of  Reference  1  except  that: 

a.  Time -dependent  quantities  used  in  Equations  [Al]  through 
[A5]  are  now  determined  by  the  digital  computer  in  accordance  with  the 
methods  of  Appendix  B  rather  than  from  data  obtained  manually  from 
analysis  of  an  oscillographic  record  of  ship  motion,  as  described  on  page 
68  of  Reference  1 . 

b.  Certain  time -independent  quantities  (i.  e. ,  hull  parameters) 

used  in  Equations  [Al]  through  [A5]  may  now  be  calculated  by  a  digital 

.  2 
computer . 
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APPENDIX  B 


DATA  REQUIREMENTS  AND  METHOD  FOR  COMPUTING  DATA 

To  solve  the  set  of  equations  ([ A1 J  -  [A5])  given  in  Appendix  A,  the 
following  time-independent  and  time -dependent  data  are  to  be  supplied  (as 
input  to  the  computer): 

Bl.  TIME -INDEPENDENT  QUANTITIES 

2 

The  hull  parameters  given  here  may  be  computed  either  digitally 

g 

or  manually  from  ship  plans.  The  time -independent  data  furnished  are: 


m  (x) 
s 

g- 

El  (x) 

h. 

KAG  (x) 

1  „  M 

HZ 

p  _  c  (x,  t) 

H  (x,  t)  cj 

Time  intervals 

The  following  quantities  required 
for  the  computation  of  P^  by 

means  of  a  digital  computer:** 


f.  Initial  conditions  for  rigid - 
body  and  flexural  motions 
4>  (0),  <j>( 0),  yh(0),  yh  (0), 

yp  (X,  0),  y(  ( X ,  0),  yf  (x,  0), 

V  (x,  0),  y  (x,  0). 

(Note:  $  is  given  in  radians) 

*  This  section  supersedes  Section  F3  of  Appendix  F,  Reference  1. 


**  C  or  k. required  for  computing  Pe  by  the  method  of  Appendix  F3, 
Reference  1,  need  not  be  provided  here.  Here  the  program  assumes  that 
the  heave  meter  is  placed  at  the  center  of  gravity  of  the  ship,  and  the 
computation  makes  use  of  that  fact.  Moreover,  (y,  z),  the  coordinates  of 
the  ship  profile  for  each  section,  are  now  required  in  place  of  c  (x),  the 
half -breadth  required  in  Appendix  F3,  Reference  1.  This  means  that  the 
points  on  the  profile  can  be  given  at  irregular  spacings,  thus  allowing  for 
a  better  description  of  the  ship  line  (discussed  later). 
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p>  U’  @s’  (y>z) 


and  the  number  of  sections  or  stations  and  their  respective  spacings  to  be 
used  in  solving  the  finite -difference  equations. 

« 

In  solving  for  (see  Appendix  B2. 5),  the  initial  values  'P  ( 0 )  and 
y^  (0)  are  required  (see  Item  f  above). 

B2 .  TIME -DEPENDENT  QUANTITIES 

Solution  of  the  finite -difference  equations*  by  means  of  the  origi¬ 
nal  AML  Code  (TRC-4)  required  that  time -dependent  quantities  (  e.g., 

77i,  *— — «}  P  ,  c,  etc.  )  be  supplied  at  several  time  intervals.  However,  an 

d  t  e 

extension  of  TRC-4  has  now  been  devised  whereby  the  time -dependent 

quantities  are  computed  vvithin  each  time  step.  This  routine  requires  only 

time -independent  data,  some  additional  constants,  and  several  points  to 

3 

describe  each  ship  line.  The  method  for  representing  a  ship  line  and  the 
method  for  determining  the  time -dependent  quantities  are  now  described. 
B2.1.  Ship  Lines 

Initially  the  revised  program  fits  a  series  of  third -degree  poly¬ 
nomials  to  the  set  of  points  describing  each  ship  cross  section  in  order  to 

3 

get  a  function  approximating  the  ship  line.  A  segment  of  a  cubic  is  fitted 
between  every  two  adjacent  points  so  that  in  a  neighborhood  of  each  point, 
the  two  adjacent  segments  coincide  (i.e.,  at  the  junction  point,  the  func¬ 
tion  and  its  first  and  second  derivatives  are  continuous).  The  resulting 
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function  may  be  written  explicitly  as  follows: 


_  _  2  o  o  .? 

f(y)  =  a  +  by  +  c  y  +  A  Q  (y-yQ  y  +  A  j  (y-y^  +  An  (y -  >2  )+  + 

•  •  •  +An-i(y-rn-i)l 


where 


y  =  0 
J  o 


is  the  first  water  level, 


y  ,  y^>  •  •  •  yM_!  are  the  second  to  N-l  water  levels  used  as 


ordinates  in  the  data  points, 
is  the  corresponding  half-breadth  at  level  y, 


(y-y^l  =  o 


=  (y-y^ 


for  y  y. 


1  =  1,2 _ y-i 


Thus  all  terms  that  include  the  factor  (y  -  y  )  are  to  be  disregarded 
whenever  that  factor  is  less  than  or  equal  to  zero. 


integrated  for  the  general  case: 


A  =J'yf(y)dy  =  (a-'£  A.y.  )  y  +  (b  +-3^/1. y{-  2)~  +  (c- 3  ^  A.y.)2_ 


+  (An  +  2  A.)l —  ,  J_  53/1.  y.4 
0  i=l  1  4  14  i=011 


*  The  notation  in  Reference  3  has  been  altered  in  the  following  discussion 
in  order  to  remain  consistent  with  the  coordinate  system  and  notation 
used  in  Reference  1.  However,  the  notation  was  not  changed  in  the  actual 
FORTRAN  coding. 
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7  K  , 

where  the  last  term  _ V  -1  •  (0  -  >•  J  =  o  for  v  <  0  .and  where 

4  jL-i  1  'i  +  Ji  = 

7  =  0 

K  =  N  -  1  when  all  the  (y  -  yj  terms  are  positive,  and  K  =  N-  J(J>  1) 

when  all  terms  (y  -  yN_  J+  j)  t0  (y  -  yN_j)  are  <  0. 

In  accordance  with  the  methods  of  Reference  3,  the  coefficients 
a,  b,  c,  A q ,  A  ,  A^  are  determined  and  stored  in  the  machine 

ready  for  use  whenever  the  function  is  needed.  In  addition,  the  number 
of  inflection  points  in  each  ship  line  and  their  locations  are  determined 
for  use  later  to  obtain  the  section  area  coefficient  17  for  the  added  mass. 

It  should  be  noted  here  that  the  function  fits  the  data  points  very 
closely,  and  care  should  be  taken  to  avoid  including  unnecessary  points 
as  input;  often,  especially  in  the  case  of  sections  near  the  midship,  in¬ 
flection  points  appear  when,  in  reality,  there  are  none.  This  generally 
occurs  when  more  than  two  points  are  used  to  describe  a  nearly  straight 
portion  of  the  ship  line.  Extreme  caution  should  be  exercised  in  approxi¬ 
mating  the  stern  section  if  it  has  a  very  narrow  portion  below  the  water¬ 
line  because  the  approximation  is  critical  in  that  region  and  the  function 
may  cross  the  axis  or  introduce  extraneous  humps. 

If  data  points  are  taken  from  the  final  ship  specifications  given  by 
the  naval  architect,  the  above  function  should  be  a  good  approximation; 
however,  if  the  points  are  not  accurate,  fewer  points  should  be  used. 
These  can  be  selected  from  a  plot  of  all  the  given  points. 
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B2.2  Determination  of  Half -Breadths  c,  c',  c"  and  Corre¬ 
sponding  Immersions  y\  y,  y1"  and  Deadrise  Angle 

When  the  dynamic  problem  begins,  the  immersion  y  is  computed 

(see  Equation  [3]  of  Reference  1): 


y  =  1)  +y  -v,  -  y 
1  w  h  p 


where  (See  Appendixes  D,  E,  and  F  of  Reference  1) 
h 


Vw  = -  COS  0„ 


ft.  - 


(V 


w 


+  U  cos  8$)  t  +  £n  cos  0g  J 


u> 


yp=(x-xc.g.W 


(v-l) 


•  (v-l\  At' 

+  r*  )-T 


<P 


(v-l) 


+ 


:m  ,  ;(v-i) 


+  P 


) 


At' 

2 


and  the  superscript  in  terms  of  v  refers  to  time,  v  being  the  current 
time  step  described  on  page  75  of  Reference  1;  although  all  the  variables 
are  functions  of  v,  only  those  equations  requiring  insertion  of  v{ov  v  ±  k, 
k  =  1,2  . . . ),  for  the  sake  of  clarity,  contain  the  superscript  here,  where¬ 
as  variables  in  the  other  equations  will  be  tacitly  understood  to  contain 
similar  superscripts  also. 

a* 

See  Notation  for  definition  of  symbols  used  here. 
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Next  the  downward  vertical  velocity  is  found  to  determine  the  di¬ 


rection  of  immersion:* 
y  =  Yw  -  yh  -  Yp 


where  ,,  9 

•  n  2n  ,,  .  „ 

y  - - V  sm  6 

w  0  ,  w  n 


h  2  n 

2  X 


V  sin 
tv 


2  nx 

n 


X 


h  2  n 

2  X 


V  sin 
w 


—  \(V 

XL1" 


+  U  cos  0  )  t  +  £  cos 
s  n 


el 

sj 


)  <A 


Y  (v)  _  9  V  — 

A  -  A 


For  each  section  the  half -breadth  c,  corresponding  to  the  com¬ 
puted  immersion  y,  is  found  simply  by  solving  the  function  approximating 
the  ship  line  at  that  immersion.  In  addition,  calculation  of  the  first  de¬ 
rivative  yields  the  slope  yS  of  the  line  at  the  point  in  question  (see 
Appendix  B  of  Reference  1). 

Szebehely's  relation,  which  accounts  for  the  piled -up  water  about 
a  body  penetrating  a  water  surface,  is  then  employed  to  determine  the 
actual  half -breadth  for  the  computation.  If  the  ship  is  immersing  (y  >  0), 
a  value  c  greater  than  the  half -breadth  c'  corresponding  to  the  ship  line 

position  will  be  determined;  conversely,  if  the  ship  is  emerging  (y  <  0), 

* 

In  the  first  computation,  the  value  for  and  y^  are  extrapolated  and 

then  recomputed  by  other  relationships  for  the  variables  (see  Appendix 
B2.5). 
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a  value  cM  smaller  than  the  half -breadth  c'  will  be  determined;  and  if 
the  immersion  has  reached  a  maximum  (or  a  minimum)  (y  =  0),  the  half¬ 
breadth  will  remain  unchanged. 

The  relation  between  half -breadths  given  by  Szebehely  is: 


where  c"<  c'  <  c,  c'  being  the  original  half -breadth,  and  tan  p  is  the 
slope  of  the  ship  line  at  the  point  (y,  c‘). 

If  a  new  half -breadth  (c  or  c" )  is  determined  from  c'  and  the 
previous  equation,  it  is  necessary  to  find  the  corresponding  draft  y'  or 
ym,  respectively.  This  is  done  with  an  iteration  based  jn  linear  interpo¬ 
lation  of  the  ship  line  function.  For  this  reason  either  the  function  must 
be  monotonic  and  single -valued  (i.e.,  corresponding  to  a  nonbulbous 
section)  or  additional  information  must  be  furnished  to  permit  determina¬ 
tion  of  which  of  the  two  possible  values  of  f  (y)  is  pertinent. 

B2.3.  Determination  of  Submerged  Area  A,  Section  Area 

Coefficient  y  ,  Added  Mass  Coefficient  Cy  and  Added 
(or  Virtual)  Mass  m 

Having  found  the  depth  of  submersion  y'  or  y'"  of  the  ship  section, 
it  is  possible  to  calculate  the  submerged  area  by  the  method  described  in 
Appendix  B2. 1  and  the  section  area  coefficient  y  making  use  of  the  inte¬ 
grated  ship  line  function  by  the  method  now  described.  The  coefficient  rj  , 
which  is  used  in  determining  the  added  mass,  is  based  on  a  modified  section 
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area.  If  the  section  has  any  concave  portions,  these  must  be  discounted 

by  drawing  a  tangent  line  between  convex  portions  of  the  ship  line  and 

7  13 

treating  the  modified  section  as  if  it  were  solid.  ’  This  is  done  on  the 
machine  by  an  iterative  process,  starting  at  the  inflection  points  (see 
Figure  1).  Holding  one  point  A  fixed  at  the  point  of  inflection  and  moving 
another  point  B  away  from  it,  the  change  in  the  slope  of  the  line  connecting 
the  points  AB  is  checked  for  a  change  in  direction  or  sign.  If  the  slope  is 
increasing  and  suddenly  begins  to  decrease,  the  point  of  tangency  (i.e., 
the  moving  point  B)  is  fixed,  and  then  point  A  is  moved  away  from  the 
inflection  point  until,  again,  the  change  in  slope  changes  sign.  This  point 
A  is  fixed  at  the  second  tangency  point.  However,  the  new  line  is  no 
longer  tangent  to  the  first  tangency  point.  Therefore,  this  process  is  re¬ 
peated  until  a  line  sufficiently  tangent  to  both  ends  of  the  concave  portion 
is  determined. 

If  there  is  only  one  inflection  point  in  the  submerged  portion  of  the 
ship  section,  a  point  on  the  upper  end  of  the  ship  line  (i.e.,  point  B)is  fixed 
at  the  waterline  and  only  the  lower  point  A  is  varied  until  tangency' is 
obtained.  If  there  are  more  than  two  inflection  points,  one  or  more  tan¬ 
gent  lines  will  be  determined  in  the  manner  just  described. 

The  area  A  is  then  found  by  adding  the  area(s)  of  the  trapezoid(s) 
bounded  by  the  tangent  line(s)  to  the  actual  area  of  the  remainder  of  the 
submerged  portion  of  the  ship  section  (see  Figure  2).  77  is  equal  to  this 
area  A  which  is  the  area  of  one -half  the  total  modified  cross  section 
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Figure  B1  -  Method  for  Bridging  Hollow  Portions  of  an  Immersed 
Ship  Section  with  Tangent  Lines 

First  A  is  fixed  at  the  point  of  inflection  and  B  moves  until  the  point  of 
tangency  pf  line  AB  to  the  upper  convex  portion  is  found.  Then  B  is  fixed 
at  this  point  and  A  moves  until  the  point  of  tangency  of  line  BA  to  the  lower 
convex  portion  is  found.  The  process  is  then  repeated  until  line  AB  is  suffi¬ 
ciently  tangent  to  both  ends  of  the  concave  portion,  i.e.,  the  upper  and  lower 
convex  portions.  Lines  AB'  and  B/l^-show  the  changes  in  slope  of  lines  AB 
and  BA,  respectively. 
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Figure  B2  -  Submerged  Area  for  a  Ship  Section  which  Includes 
a  Hollow  Portion 
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divided  by  the  product  of  the  half  breadth  and  the  draft  at  the  water  level 
(i.  e. ,  c  •  y'  or  c"  •  y,M  where  c,  c"  and  y\  y"1  are  determined  in  the 
manner  just  described);  i.e., 

A  A 

V  =  - ;  or 

c.y  c  y 

The  added  mass  is  now  determined  in  the  manner  described  by 
7  1  ?  4  5 

Lewis,  Prohaska,  ’  and  Landweber.  ’  For  a  cylindrical  hull 
partially  submerged  and  vibrating  vertically,  the  added  mass  per  unit 
length  is  given  by 


n  2  r 

m  =~2~P  C 

where  p  is  the  density  of  sea  water, 

c  is  the  local  half -breadth  at  the  waterline;  in  general,  c  -  c 
or  c"  here  (see  Appendixes  B  and  F3  of  Reference  1),  and 

C^is  a  coefficient  depending  upon  the  form  of  the  submerged  cross 
section  determinable  as  a  function  of  17  ,  c,  and  y'  (or  77,  c", 
and  y'"  if  the  ship  is  emerging  rather  than  immersing). 

The  inertial  coefficient  Cv  may  be  determined  by  a  conformal 

transformation  of  the  known  flow  about  an  elliptic  cylinder  of  approximate 

ship  section.  Consider  the  mathematical  transformation  used  by  Land- 

4  5 

weber  and  Macagno  ’ 


Z  =  .  +^i.+-5L+iL 

zm  zn  zp 


ar  a3 ’  °5  real 
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in  which  the  indices  m,  n,  etc.,  are  odd  numbers  for  the  case  of  symmet¬ 
rical  sections.  This  transformation  comprehends  the  transformations  with 

indices  (m,  n)  =  (1,  3);  a^  =  0,  and  (m,  n)  =  (1,5)  (1,7)  (3,  7)-  a^  =  0  used  by 
6  7  13 

Lewis  and  Prohaska,  ’  respectively,  to  develop  an  analytical  pro¬ 
cedure  for  evaluating  the  virtual  mass  of  a  class  of  two-dimensional  forms, 
representative  of  ship  sections,  by  means  of  conformal  mapping  of  a 
circle.  Here  z  is  the  complex  variable  of  the  circle  plane  and  Z  is  the 
complex  variable  of  the  transformed  plane.  Different  values  of  a^,  a^, 
and  a  give  different  sectional  forms  and  corresponding  variations  in  the 

D 

two-dimensional  virtual  mass. 

For  the  three-parameter  family  of  two  -dimensional  forms  (a,_=£  0) 
with  (m,  n,  p)  =  (1,  3,  5),  Landweber  gives  the  following  results  (see 
Equations  [17],  [16]  (alternate  form  above  it),  [4],  and  [5]  of  Reference 


D  =' 


y / 
$ 


1  +  Qj  +  Ctg  + 

l-aI+a3-a5 


A  „  (1  -a2j  -3  a3  -5  a5)  (1  -  a*  •  3  <*3  •  5a  g) 

2cy'  4  (1 +ci][+a3+a5)(l-a1+a3-a5)  4  (J  +  j2  _  (a  ^  +  a  _  ;  2 

(Only  immersion  variables  c,  y'  are  considered  here.  Emersion  vari¬ 
ables  c",  y'"  are  treated  similarly.) 
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Special  cases  of  these  results  are: 

Lewis  Forms:  (m,  n)  =  (1,  3);  a  =  0 

o 

Prohaska  Forms:  (m,  n)  =  (1,5),  (1,7),  (3,7);  a  =  0 

o 

For  the  Lewis  forms,  a^  and  a^  are  solved  in  terms  of  p  and  17, 

which  are  known  (predetermined)  quantities  from  the  last  two  equations 

(a.  =  0).  Thus 
0 

*1  =*■ 

4  (p  n  +  [>n-T7T  +  4l>ri ) 

1  -  p  \(p+l)  a 

a  =  _ 

3  p-1 

J 

Since  a  and  a  each  assume  two  values,  it  is  necessary  to 

-L  o 

establish  some  restrictions  on  the  parameters.  Consider  that  only  the 
function  z  =  x  +  iy,  which  lies  in  the  first  quadrant,  yields  useful  forms 
lying  within  the  rectangle  of  width  2c  and  depth  y'.  For  this  case,  the 
following  restrictions  given  by  Landweber  must  be  imposed  (see  Equation 
[41]  of  Reference  4): 

H  1  y' 

*  Landweber 's  \  = - ►  —  =  _ 

B  p  c 
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1 

iT*  ’ 


1 


L  _L  <  J_  <  J_  j_  ._L  > 

l>  2  c  [>  +  -<  '  1> 

Only  one  set  of  values  for  and  a^  meet  these  requirements, 

«•* 

hence  a  unique  value  for  C^.  can  be  determined. 

The  coefficients  for  the  Prohaska  forms,  a.  and  a_  or  a,  and  a_ 

1  5  17 

or  a  and  a  ,  can  be  solved  in  a  similar  manner.  These  forms  as  well 

O  i 

as  the  Lewis  forms  are,  however,  a  special  case  of  the  Landweber  forms 
now  to  be  treated. 

To  solve  the  Landweber  three-parameter  family  of  forms  for  a^, 

a^,  and  a^,  we  first  introduce  the  parameter^ 

I 

1  "  - 3 

c(y'/ 

where  I  is  the  moment  of  inertia  of  the  submerged  cross  section  about 
the  transverse  axis  in  the  free  surface 

y  =  y'  2 

l  =  2  f  f  (y)  y  dy 

I  and  I  are  known  (predetermined)  quantities. 

*  The  resulting  added  mass  per  unit  length  m  must  be  divided  by  2000  lb 

2  2  2 
and  32  ft/sec  in  order  to  have  m  in  terms  of  ton-sec  /ft  ,  the  unit  of 

mass  in  Reference  1. 
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Now  a^,  a^,  and  are  solved  in  terms  of  p,  77  ,  and  I  (all  known 
quantities)  using  the  equations  for  p  and  n  given  previously  and  the  follow¬ 
ing  equation  (see  Equation  [10]  of  Reference  5);  the  details  of  solution  are 
not  given  here.  The  results,  of  course,  include  those  of  Lewis  and 
Prohaska  as  special  cases. 

-  2  22 

ao  +  2/3  (a  -  a  a  +  a  )  -  a  a  ] 
3  3  35  5  35 

-  2 

)  +12/3  a 


1  = 


'(f)4 


|  2  [f  +  f 


1  28 


-4  -3-2  2  2 

-  B  -2/3  a  -  2R  (4a  -5a  a  +6a 

3  3  3  5  5 


2  2  2  2 
-  (3  a  +  a  )  (a?  +  0  a 

356  0 


here  B  =  1  -  a  . 

1- 

To  meet  the  practical  requirement  that  the  forms  should  be  within 


a  rectangle  of  width  2c  and  depth  y',  the  following  conditions  and  restric- 
tions  are  imposed: 

2 

Conditions  for  ct  =  —  • 


Lower  limit: 


r  >  r 
1  <  2 


for  a_  <  £ 

O  ^  j 


Upper  limit: 


S1  <  S2 


for  a  5  — 

5  ^  27 


*  Calculations  of  the  added  mass  using  these  relations  have  been  made  on 
the  IBM  digital  computer  at  the  Computer  Center,  University  of  Califor- 

14 

nia,  Berkeley,  California.  The  program  for  the  control  of  the  auto¬ 
matic  computations  permits  selection  of  any  foreseeable  number  of 
coefficients — many  simple  calculations  have  been  done  with  various  com¬ 
binations  of  coefficients  through  a^ . 

**  °^and  o  (page  are  notations  in  Reference  5  equivalent  to  notations 
2/c  and  n  respectively,  in  present  report. 
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To  ascertain  which  forms  associated  with  a  particular  class  of 
ship  (or  even  sections  of  a  given  ship)  are  preferable  —  as  determined 
by  comparison  between  the  theoretical  and  experimental  response  —  it  is 
desirable  that  C^.  be  computed  for  the  Lewis,  Prohaska,  and  Landweber 
forms;  the  latter  two  for  various  pertinent  confirmations  of  m,  n;  and 
m,  n,  p,  respectively.  This  would,  of  course,  only  require  coding  of  the 
Landweber  forms  inasmuch  as  the  other  fcrms  are  included  as  special 
cases.  Reference  1  used  Figure  24  of  Reference  7,  corresponding  to 
(m,  n)  =  (1,  5),  to  compute  C^..  The  final  results  for  the  time-varying 
stress  on  the  hull  girder  amidships  using  values  of  so  computed  were 
in  good  agreement  with  experiment  (see  Figure  23  of  Reference  1).  Good 
judgment  in  the  initial  choice  of  the  analytical  form  to  be  used  for  a  par¬ 
ticular  ship  cross  section  (which  may  vary  between  the  extreme  U  to 
extreme  V  in  actual  form)  may  be  applied  by  comparing  the  theoretical 
values  and  the  experimentally  obtained  values  of  Cy  or  virtual  mass  for 
the  forms  (U,  UV,  V)  and  indices  (m,n)  used  in  References  7,  13,  and  15. 


dx 

B2.4.  Determination  of  Velocity  of  Advance  and  Downward 
Relative  Velocity  y 


The  downward  relative  velocity  of  the  ship,  y  ,  is  now  found 
(see  Equations  [2 j  and  [3]  and  Appendixes  D  and  E,  Reference  1,  and 
Notation): 

yr  =  yw  '  *h  "*P  +  (u  -u  cos  V  * 


dx 

=  y  "  dF 


0 
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where 


(U  -u  cos  6  ) 
s 


|  ^JL  v  cos  e 
2  A  w  n 


B2.  5.  Determination  of  Forcing  Function  P  and  Rigid -Body 

•  ^ 


Motions  yh,  yh,  y 


P  =  P  -  gm  =  -j—  (my  )  +  (g  +  v)  P A  -  gm 
ee  sdt  r  s 


where 


=  y  =  -iL  (llL\ 

yw  2  y  A ) 


V  (V  +  U  cos  0  )  cos  6 

w  w  s  n 


and  where 


9  ” 

,  =  fJL  (v  +  U  cos  0  )  t  +  <f  cos  0 

n  a  w  s  n  s 


The  first  computation  of  P^  may  not  be  very  accurate  due  to  its 
dependence  on  extrapolated  values  for  *P  and  y  .  Hence,  the  program  re- 
computes  values  for  'P  and  y^  as  follows:  (See  pages  68  and  76  of 
Reference  1.) 
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/_  \  (v-%) 

lo  (  ">  A-)  .»  ; 

'V  (. .,) 

n=0  \  s  J  n+'/2 


■  M-  ■  ("-'I  + 

h  '  \ 


.  (v)  .  (v-1) 

</'  =  <!> 


N:1  /_  \  -  ’A) 

»=o  \  e  J  k  +  >a  (Xn+%‘Xc-S‘) 

,1°  ("si’)»  +  K  (*-+»■’'«) 


where  x  is  the  distance  from  the  center  of  gravity  to  the  stern  and 
c.g. 

the  subscript  n  +  1/2  refers  to  the  midstation  values  between  Stations 
n  and  n  +  1 . 

Other  parameters  are  computed  from  these  new  values: 

(v)  (v-1)  .(v)  .(v-1)  At' 

h  ~yh  +(yh  +  yh  }~r 


t(vKt(v-l)+^vKi(v’l)y)^f 


y(p)  (*  ,i/~x  ) 

J  t>  r  » +  72  C.g. 


(v)  (v-1)  (v)  (v  1)  A  / 

y„  =y  +  <y„  +  y„  ) 

p  p  P  P  2 
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Subscripts  on  y  ,  y  are  tacitly  understood  to  pertain  in  the  actual 
P  P 

coding.  Then  the  immersion  y1  is  recomputed  from  the  newer  values  of 
its  dependent  variables  and  the  subsequent  computations  are  repeated. 
The  resulting  force  Pg  is  compared  to  the  first  Pg.  If  they  are  not  suffi¬ 
ciently  close,  then  y^  and  are  again  recomputed  and  the  process  is 
repeated  until  a  satisfactory  P^  is  found  for  each  section. 


B2.6.  Determination  of  Total  Mass  per  Unit  Length  /*  (x,  t) 
{Including  an  Allowance  for  Virtual  Mass  m  (x,  t)  and 
Damping  Force  <?(x,  t)) 

dx  — 

Thus  the  time -dependent  parameters  m,  and  P^  are  used  to 

determine  u(-m  +  m  )  and  other  variables  (e.g.,  P)  required  to  find  the 

s 

transient  vibration  of  a  ship  using  the  equations  in  Appendix  A.  In  par¬ 
ticular,  c  ,  the  time -dependent  damping  force  per  unit  velocity  per  unit 
length,  is  found,  from  the  expression  .5.  =  K.  The  constant  K  and 

[t  ( Xy  t  )Ct) 


modal  frequency  w  are  given  (cjcan  be  determined  easily  by  obtaining  a 

g 

normal  mode  solution  using  a  code  devised  by  the  Applied  Mathematics 
Laboratory  of  the  Model  Basin  or,  more  roughly,  from  Schlicks  Formula). 


The  value  of  the  constant  K  depends  on  the  class  of  ship  and  is  determined 
in  accordance  with  the  methods  of  Reference  16;  w  =  1  if  it  is  desired  to 
use  the  damping  constant  as  a  frequency  independent  parameter. 

The  remainder  of  the  program  is  essentially  unchanged  from  the 
original  TRC-4  used  in  Reference  1. 
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APPENDIX  C 


COMPARISON  OF  PRESENT  METHOD  OF  COMPUTATION 
WITH  ALTERNATIVE  METHOD  GIVEN  IN  APPENDIX  F3  (PAGES 

80  -  82)  OF  REFERENCE  1 

It  is  of  interest  to  compare  the  present  method  with  the  alterna¬ 
tive  method  of  Reference  1  (pages  80-82). 

The  present  method  uses  the  coordinates  to  determine  a  function 
approximating  the  ship  profile  (i.  e. ,  a  function  of  half -breadth  versus 
immersion)  thereby  superseding  the  procedure  described  on  pages  80-81 
of  Reference  1  (paragraph  a)  for  determining  the  half-breadth  for  a  given 
immersion  using  a  data  storage  technique.  Only  one -half  the  ship  pro¬ 
file  need  be  considered,  due  to  symmetry. 

The  present  method  uses  the  ship  line  function  in  the  dynamic 
problem  to  find  a  value  for  fi  by  taking  the  arctangent  of  the  first  deriv¬ 
ative.  This  supersedes  the  procedure  described  on  pages  81-82  of 
Reference  1  (paragraph  b)  for  determining  by  interpolation. 

The  present  method  obtains  the  submerged  area  A  by  integration 
rather  than  by  using  an  approximating  sum  as  described  on  page  82  of 
Reference  1  (paragraph  c). 

In  accordance  with  the  present  method,  the  added  mass  is  com¬ 
puted  at  the  time  and  immersion  when  it  is  needed  and  not  prior  to  the 
time -dependent  calculations  as  described  on  page  82  of  Reference  1 
(paragraphs  d  and  e). 
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The  computation  of  the  immersion  and  the  hydroaynamic  force, 
etc. ,  is  described  in  Appendix  B  of  this  report. 
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APPENDIX  D 


EXTENSION  OF  ANALYSIS  TO  IRREGULAR  WAVES 


The  wave  profile  treated  in  this  analysis  is  a  regular  wave;  see 

page  51  and  Appendix  D  of  Reference  1.  Such  simple  waves  can,  however, 

17 

serve  as  a  basis  for  the  synthesis  of  more  complex  waves.  For  ex¬ 
ample,  an  irregular  wave  moving  in  one  direction,  with  wave  crests 
parallel,  may  be  treated  as  a  linear  summation  of  regular  sine  waves  of 
different  wave  lengths  and  periods.  This  representation  can  be  extended 
by  adding  a  different  constant  phase  angle  to  each  component  wave.  Thus, 
the  treatment  given  in  this  report  for  a  regular  wave  can  be  extended  tc 
include  the  terms  necessary  for  a  deterministic  representation  of  an  ir¬ 
regular  wave  (i.e. ,  arbitrary  wave  of  encounter  for  a  unidirectional  sea). 

An  irregular  wave  can  also  be  characterized  from  a  statistical 
point  of  view.  Thus  an  element  of  probability  may  be  introduced  into  the 
description  of  the  sea  by  statistically  distributing  the  phase  angles  in  ac¬ 
cordance  with  a  certain  distribution  function  so  that  all  phase  angles  are 

17  17 

equally  probable.  For  this  case  the  spectral  density  function  repre¬ 
sents  the  manner  in  which  the  energy  of  the  composite  wave  is  distributed 
over  the  frequency,  thereby  fully  characterizing  the  nondeterministic 
irregular  wave  of  this  type. 

The  sea  waves  may,  however,  exhibit  an  irregularity  in  two 
directions.  These  irregular  waves  may  be  synthesized  from  component 
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sine  waves  assumed  to  move  in  any  direction.  The  spectral  density 
function  now  represents  the  manner  in  which  the  energy  of  the  composite 
waves  is  distributed  in  both  the  frequency  interval  and  the  direction  in¬ 
terval.  Extension  of  the  treatment  given  in  this  report  for  a  regular  wave 
may  yield  a  practical  procedure  for  treating  nondeterministic  (statistical 
or  probability)  representations  of  irregular  waves. 
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